Non-equilibrium Characterization of Spinodal Points using Short Time Dynamics 
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Though intuitively appealing, the concept of spinodal is rigorously defined only in systems with 
infinite range interactions (mean field systems). In short-range systems, a pseudo-spinodal can be 
defined by extrapolation of metastable measurements, but the point itself is not reachable because it 
lies beyond the metastability limit. In this work we show that a sensible definition of spinodal points 
can be obtained through the short time dynamical behavior of the system deep inside the metastable 
phase, by looking for a point where the system shows critical behavior. We show that spinodal points 
obtained by this method agree both with the thermodynamical spinodal point in mean field systems 
and with the pseudo-spinodal point obtained by extrapolation of meta-equilibrium behavior in short 
range systems. With this definition, a practical determination can be achieved without regard for 
equilibration issues. 
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Keywords: spinodal; short time dynamics; Monte Carlo; Potts and Ising models 



I. INTRODUCTION 

First-order phase transitions are accompanied by hys- 
teresis and metastabiUty: even though the thermody- 
namic transition happens at the value (pt of the control 
variable ip, when (p is varied smoothly from ipi > Lpt to 
Iff < ipt, the system remains in the phase correspond- 
ing to thermodynamic equihbrium a.t ip > (pt (and con- 
versely when changing ip in the opposite sense). When 
the phase survives being carried out beyond its thermo- 
dynamic "homeland" , it is called metastable. Metastable 
phases have a finite lifetime, but this time can be very 
long. Diamond at room temperature and pressure, and 
glass- forming supercooled liquids are very well-known ex- 
amples of long-lived metastable phases (so long-lived, in 
fact, that for many purposes they can be considered an 
equilibrium phase). In general, however, the metastable 
phase cannot exist for all and it is not observed if v/ is 
less than some value ips-p- This is the idea behind the con- 
cept of spinodal point. However, to define the spinodal, 
some care is required. 

At the mean-field level, the spinodal is well de- 
fined. Focusing on the ferromagnetic case to be specific, 
let's consider the extended free energy per particle [ij 
/3(r, m, /i), with h the magnetic field and m the magneti- 
zation, /a depends on two conjugate variables because it 
is defined such that the probability of finding a value M 
of the magnetic moment is cx exp[— /3A^/3(r, Af/A^, ft,)], 
where N is the system size. In mean field, and in the 
limit N oo, fs has two minima as a function of m 
for T below some critical ternperature Tc and h within 
some range —hgp < h < hgp P, [I]- As a function of h, 
the (first-order) transition occurs ai h = he = 0. When 
h — 0, the two minima are symmetric, corresponding to 
the two broken-symmetry phases. When ft 7^ 0, the abso- 



lute minimum corresponds to the thermodynamic equi- 
librium (or stable) phase, while the secondary local min- 
imum defines a phase which is dynamically stable but 
of higher free energy: the metastable phase. When it 
exists, the (mean field) metastable phase has infinite life- 
time. At ft = ftgp, the secondary minimum disappears 
(it becomes an inflection point) and the phase becomes 
unstable: a system prepared with a magnetization differ- 
ent from the (thermodynamic) equilibrium value evolves 
toward the equilibrium state. For |ft| > ftgp, /a has only 
one minimum. Thus in mean field the spinodal, which is 
the point where the metastable phase becomes unstable 
(in the sense that a susceptibility becomes negative), is 
also the limit of metastability, i.e. the point up to which 
the metastable phase can be observed. 

When the interactions have a finite range, matters are 
more complicated 0, Q. On one hand, the metastable 
phase ceases to be observable before it becomes unsta- 
ble ^, ^] . This is because as the system moves away 
from the transition, the lifetime of the metastable phase 
decreases while its relaxation time increases. When they 
become of the same order, the phase is unobservable. 
This is the metastability limit, which is thus different 
from the spinodal. The metastability limit is also called 
kinetic spinodal Q, and the term thermodynamic spin- 
odal is sometimes used for the spinodal as defined above 
(onset of an instability). On the other hand, the order 
parameter can fluctuate in space. Although /s can still 
be deflned (and can be used to compute the true equilib- 
rium properties), in the thermodynamic limit it has no 
convexity changes and only one minimum; therefore, it is 
useless to define a spinodal. Further, although the exten- 
sive F3 has a double- well shape [8], the local maximum 
cannot be interpreted as a barrier to the growth of the 
stable phase, because the global magnetization is not a 
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good coordinate to describe this process: a supercritical 
nucleus of the stable phase (one whose growth is ther- 
modynamically favored) can form without change in 
global magnetization [8]. The eventual disappearance of 
the secondary minimum in is hence unrelated to the 
spinodal. 

The spinodal is then beyond the metastability limit, 
and hence outside the realm of thermodynamics and of 
quasi-equilibrium treatments. Due to these difficulties, it 
has been concluded that the spinodal only makes sense 
in mean field [lOi] or in finite size systems ^3;]. However, 
signs of an instability are detectable in (meta)equilibrium 
measurements: the susceptibility and relaxation times of 
the metastable phase increase as one goes deeper into 
the metastable region, and if extrapolated with a power 
law, seem to diverge at a point beyond the metastability 
limit, called pseudospinodal [loL 

Here, rather than to abandon the idea of a spinodal 
in finite-range systems, we propose to define it through 
out-of-equilibrium properties. The idea is that the spin- 
odal should be a point with infinite susceptibility and 
infinite relaxation time. In this sense it resembles the 
critical point of a second-order phase transition. It is 
known fl3, ^T§\ that it is possible to detect critical points 
by studying the short-time dynamic behavior of the or- 
der parameter and correlation functions. We propose to 
use the same method to identify a point deep into the 
metastable region which, at short times, behaves dynam- 
ically like a critical point. This pseudocritical regime (in 
the sense that it lasts only for a finite time) does not im- 
ply the existence of a thermodynamic singularity (see dis- 
cussion in sec.|V]). We show that in a mean-field model, 
where the spinodal is well-defined and can be worked out 
analytically, the point identified with this technique is 
precisely the thermodynamic spinodal. Thus by defining 
the spinodal as the point where this pseudocritical dy- 
namics takes place we provide a sensible generalization 
of the spinodal concept to finite-range systems. We apply 
the method to the two-dimensional g-state Potts model 
[T^ with q > 4, where it gives a reasonable result, provid- 
ing a bound for the metastability limit and locating the 
spinodal very near to the pseudospinodal. The technique 
has the advantage that it does not need equilibrium data, 
which is an essential requirement in order to determine 
the pseudospinodal. 

Our proposal was inspired by the finding of Schiilke 
and Zheng [l5| that the short-time dynamics (STD) 
applied to the Potts model defines two "pseudocritical 
points", which are closer together the weaker the (first- 
order) transition, and coincide for second-order transi- 
tions. A similar situation was observed in models with 
out-of-equilibrium transitions [H, [O] ■ 

The paper is organized as follows. In Sec. [H] we briefly 
review the STD technique. In Sec. lIIII we consider mean- 
field spin models and we show that the STD method ac- 
curately describes the thermodynamic spinodal in those 
models. In Sec. II VI we apply the method to the ferromag- 
netic q-state Potts model with nearest-neighbor interac- 



tions in two dimensions. Sec. |V] presents our conclusions. 

II. USING SHORT-TIME DYNAMICS TO 
IDENTIFY CRITICAL POINTS 

The STD technique to identify critical points has been 
reviewed in Refs. 13, 221. Briefly, it is based on the 
work of Janssen et al. who studied model A (a 0'* 
Hamiltonian with Langevin dynamics) in the out of equi- 
librium regime where correlation functions are still non- 
trivial functions of two times and the order parameter 
is still time-dependent. For the present work, the rel- 
evant result is that the n-th moment of the order pa- 
rameter m'^"^(t) = (^[m(t) — (m(t))] ) obeys the scaling 
formflllll] 

TO(")(i, r, L, mo) = b-^'^/'gnib-H, b^'^r, L/b, 6^mo), 

(1) 

where t is time, t is the reduced temperature t = (T — 
Tc)/Tc, L is the system size, toq is the initial value of the 
order parameter (assumed nonzero but small), and 6 is a 
rescaling parameter. ^ is a new universal exponent that 
describes the short time behavior, while /3, v, and z are 
the usual critical exponents 

From Eq. ([T]), setting b — t^^^, for large values of L 
and small values of t^^^niQ one obtains 

mit,T,mo) ^mot^Fit^/^^T), 9^^^^^, (2) 

so that precisely at the critical point t — Q, the order pa- 
rameter obeys a power law m(i) ~ . Similarly, setting 
Too = 0, one gets for the second moment at the critical 
point 

TO(2)(t) ^t-i/^-W-. (3) 

It is generally assumed (and in agreement with numer- 
ical results [13I. [2^) that when the initial condition is the 
ordered state (too = 1), a scaling similar to Eq. ^ holds: 

m(")(i,T,L) ^ b-''^'''gn{b-H,b^'''T,L/b). (4) 

From this equation one obtains, for too — 1 and large L, 

m{t) = r^/^'Git^/^'T), (5) 

and taking the derivative of log to, 

d\ogm{t,T) ^ ^1/^^ 

The critical point can then be obtained by perform- 
ing Monte Carlo simulations at several temperatures and 
looking for the value of T at which the power laws in 
time predicted by Eqs. ^ and ([5]) hold. In addition, 
these equations together with Eq. ([6]) allow to determine 
the critical exponents p^ . 

Here we apply the above method to look for singular 
behavior in the metastable region of a first order phase 
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transition. By tuning the appropriate control parameter 
(external field or temperature) we look for a value where 
the power laws ([2]), ([3]) and ^ hold at short time scales 
(for very long times, the approach to the proper equilib- 
rium phase is seen). This value of the control parameter 
can be sensibly defined as a spinodal, as we show in the 
cases studied below. 



III. SPINODAL POINTS IN SYSTEMS WITH 
LONG RANGE INTERACTIONS 

In this section we consider spin models where each spin 
interacts with every other spin. For those systems mean 
field theory is exact and therefore provides a first check 
for the STD method. Let us first briefly review the ana- 
lytical mean-field results on the spinodal points. 



A. Thermodynamic spinodals 

a. Curie- Weiss-Ising model. We first consider the 
Curie-Weiss, or fully-connected, version of the Ising 
model. In the presence of an external magnetic field h 
the Hamiltonian is given by 



n. 



cwi 



N 



(7) 



where TV is the number of spins (s^ = ±1), h is an external 
magnetic field and J > 0. The extended free energy per 
spin, /s, can be computed exactly because Hew is an 
explicit function of the total magnetization M = Si, 
and the partition function can be evaluated for fixed M. 
The result, in the limit iV — s- cx), is 
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1 + m , 1 + m 1 — m , 1 — m 
In 1 In 



-m — /irn, 



(8) 



where (3 — 1/T (we take Boltzmann's constant equal to 
1) and m — M/N. The absolute minimum of with 
respect to m defines the stable (equilibrium) solution 
m{T^h). The model shows a second order transition at 
h = he = and critical temperature Tc = J. For T < Tc, 
there is a line of first order transitions at /i = 0, where 
m{T, h) is singular. However, an analytic continuation 
m+{T, h) from positive to negative h exists as long as 
is not too big (and conversely a continuation to_(T, /i) 
from negative to positive h) . These continuations are the 
metastable states and correspond to local minima of /s. 
Thus the conditions 



dh{T,m,h) 
dm 



= 0, 



dim? 



> 0, (9) 



define the (meta)stable states. The secondary minimum 
(and thus the metastable solution) ceases to exist when 



h < hip^ 



~hsp, the spinodal field, given by 
d^f3{T,m,h) 



dm? 



^ 0. 



(10) 



Since d'^f^/dm^ — Xt^ i the susceptibihty diverges 
at hip , and it is straightforward to show that for 



h - hip^ < 1 



XT 



(11) 



In fact, this singularity can be treated like a usual critical 
point. For example, for fixed T we have 



C,^(/i-/i(-))-i/2 
Am^ih^hi-^)'/' 



where C'h is the specific heat and Am = to 



,(-) 



(12) 
(13) 

, with 



,(-) 



m{hip'') (note that Am > 0). Fixing h, an 



expansion in < Tgp — T ^ 1 gives 

Am ~ (T,p - T)i/2 
C„~(T,p-T)-i/2 

XT - (T.p - T)-i/2 



(14) 
(15) 
(16) 



If we choose Ato as order parameter, the singular be- 
havior in the neighborhood of the spinodal point can be 
characterized exactly as in a true critical point. The set 
of critical exponents is 

/5=l/2, a =1/2, 7 = 1/2, (5 = 2, (17) 

which satisfy Rushbrooke and Widom scaling relations 



a 2/3 + 7 = 2, 

7 = /3((5-l). 



(18) 
(19) 



If we assume that 77 = like in mean field critical points, 
from the Fisher scaling relation, 7 = v(2 — f]), we have 
V — 1/4. Finally, from the Josephson hiperscaling. 



vd 



a. 



(20) 



we can guess a critical dimension rfc = 6, a result con- 
firmed by a renormalization group analysis [23j . 

b. Curie-Weiss-Potts model. The Curie- Weiss- 
Potts model is defined by the Hamiltonian 



N 



(21) 



where J > 0, ct^ — 1, 2, . . . , 5 and 5{ai^ Uj) is the Kro- 
necker delta. The extended free energy per spin for 
this model can also be computed exactly |l3| 
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h{m,T) 



J(g-1)^ 
2q 



[l + iq- 1) m] In 



1 + (g — 1) m 



+ (q - 1)(1 - rn) In 



(22) 



The order parameter m is defined as 
1 



obtained from Eqs. ([9]) and pOl 



9-1 



(g (max(mfe)) - 1) . 



(23) 



N 



where = '^("'ii k ^ 1. 



and the max- 



imum means choosing the value of k that gives the 
highest value of m^. The extended free energy ([22]) 
presents two local minima (and therefore a first order 
transition) for g > 3, where the transition happens at 
Tt = J{q - 2)ln(g - l)/2{q - 1) [24]. A disordered 
metastable solution with m = (supercooled paramag- 



net) exists for rip ^ <T <Tt (with T^V ' = J/g), as well 
as an ordered metastable solution with m = mgp =/= for 

Tt < T < Tsp '' , where the spinodal temperatures are ob- 
tained from Eqs. pil| with h — 0. An analysis similar to 

the above shows that close to the upper spinodal Tip''-', 
the order parameter Am, the specific heat and the sus- 
ceptibility show the same critical behavior Eqs. (|14p - (fT6)) 
with the same critical exponents P?|) . Of course, the 
same set of exponents is obtained from other mean field 
approximations: the Landau cj)'^ model and — for the up- 
per spinodal point — the Landau cj)^ model. The values 
of these exponents are confirmed by means of renormal- 
ization group techniques 23, 25]. 

On the other hand, the behavior close to the lower 
spinodal temperature rip ■* is different. The susceptibil- 
ity shows a power-law divergence with exponent 7 1, 
while the magnetization and the specific heat remain fi- 
nite (in fact, they are identically zero in the disordered 
phase). The vanishing of the specific heat is peculiar 
to the approximation; the point is that in principle one 
should not expect to find the same critical exponent in 
both spinodals. 



(-) 



B. Short time dynamics 

We now apply the STD procedure to the Curie- Weiss- 
Ising model in the neighborhood of the spinodal point, 
using Am as an order parameter. We consider the pro- 
cess at fixed temperature, with the magnetic field as the 
control variable. We choose the initial condition in the 
ordered state (m = 1) corresponding to /i — > 00. In this 
case, the procedure analog to a quench T near Tc for a 
thermal second order transition is to set hto a value near 
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hsp ■ To determine hsp 



hsp )/^sp . 



we assume a scaling like 



h = 

Eq. H]) with T ^ {h 

In this test case, we can check the spinodal field and 
magnetization found with STD against the exact values 



"sp 



T 



f3hip^ 



■In 



1 



^sp 



T 



(24) 
(25) 



We have simulated the dynamics of this model using 
a standard Metropolis algorithm in a system of iV = 
1.6 X 10^ spins. We started with all spins up and did 
n ^ 10'^ runs while recording Am and the second moment 
of the magnetization per spin, Am'^^. Figure [T] shows 
the short time behavior of Am and Am'^^-' for T — ^Tc, 
where time is measured in Monte Carlo Steps (MCS, one 
MCS is defined as a full cycle of N spin update trials). 
From Eqs. ^ and ^ we have phip^ -0.714627 



and rri: 



(-) 

sp 



0.745356. 



For both the first and second moment (Eqs. ([3]) 
and (O) a power law behavior Am ^ t^^ and Am'-^-' ~ 

V at h ^ hip'' is found, with exponents x = 0.98 ± 0.02 
and y = 1.03 ± 0.02 respectively. The power laws are 
rather short-lived (lasting up to t ^ 100 MCS), and both 
observables deviate afterwards. However, as we will show 
below this is a finite size effect, and a true power law 
should be observerd in the thermodynamic limit. In or- 
der to calculate the derivative of log Am from simula- 
tions we have taken the magnetic field very close to the 
spinodal field, that is h — hgp ± e with e = 2 • 10^"*. Fig- 
ure [5] shows the numerical derivative obtained from runs 
at three values of h. Again these data can be fitted with 
a power law with exponent w — 2.02 ± 0.02. 

Using the scaling relations Eqs. (O and © from sec.HIl 
and the values of x and w, we obtain 

/3 = 0.49 ±0.01, (26) 

and assuming dc = 6, from Eqs. ^ and ^ we have 



z = 2.01 ±0.03. 
The exponent ly is obtained from Eq. ^ 
V = 0.249 ±0.004. 



(27) 



(28) 



Thus, in this test case we find that the spinodal field as 
well as the (static) critical exponents determined through 
STD are in excellent agreement with the theoretical re- 
sults PT ]) . (jM)) and (P5|) . 

We have also explored finite-size effects, which we use 
to obtain another determination of z. Figure [3] shows 
the evolution of Am for systems with number of spins 



ranging from iV = 1.6 x 10^ to iV = 1.6 x 10*^ at = ft.; 
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FIG. 1: (Color online) Short time behavior of the first and 
second moment of Am for the Curie- Weiss-Ising model. Full 
(black) lines correspond to power law fits at ft = —0.7146 — 



The same power law fit of Figure[T^ is also included here. 
There are clear finite size effects, and it is seen that the 
power-law fit is valid for t > tmic 10 up to a time which 
increases with system size. The time t^i^ is a microscopic 
time scale and it is the time required sweep away the 
microscopic short-wave behavior [22|. We define a time 
scale t* , so that for t > t* Am has abandoned the power 
law and is rapidly evolving toward its equilibrium value 
(here we take t* as the time when Am = 3 x 10^^). Fig. [3] 
(inset) shows t* vs. iV in a log-log scale. Such behavior 
can be understood by considering the evolution of the 
correlation length ^. In the critical regime, it increases 
following the relation 



C oc t 



l/z 



(29) 



We expect deviations from this critical law when the cor- 
relation length becomes of the order of the linear dimen- 
sion of the system, ^ ^ L. If we assume that our sys- 




10 

time [MCS] 

FIG. 2: Short time behavior of the derivative of log Am with 
respect to the reduced field t = [h — /isp '/''•sp') evaluated 
at r = 0; the straight line corresponds to a power law with 
exponent w — 2.02. 



tem of N spins behaves like one with linear dimension 
L = N^/'^", the deviation would appear when ^ ~ 7V^/'''=, 
or {t*y^^ cx N^/'^", and therefore 



t* cx TV", 



z 
dc 



(30) 



In this way we can use finite size effects to extract z. 
From the fit of Fig. [3] (inset) we get u = 0.34 ± 0.01, 
which together with dc = & gives z — 2.06 ± 0.06, in 
excellent agreement with the previous estimate. 

We conclude from this test that the STD technique 
identifies the (thermodynamic) spinodal points consis- 
tently with the static results. We proceed now to a short- 
range model, where static approaches are unsuitable. 



IV. SPINODAL POINTS IN THE 
SHORT-RANGE POTTS MODEL 

We consider now the nearest neighbor g-state Potts 
model on the square lattice, with Hamiltonian 



(31) 



<hj> 



where the sum runs now over all the pairs of nearest- 
neighbor sites. The 2-d Potts model undergoes a second 
order phase transition for q ~ 2,3,4 and a first order one 
for q > 4. In the square lattice, the transition tempera- 
ture is known exactly to be[2^ Tt{q)/J= 1/ In (l -I- ^) . 
For q larger but near q = 4, however, the transition is 
very weak. As Binder [27| pointed out, the pseudospin- 



odal temperatures Tij^-* and Tgp ' are extremely close to 
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FIG. 3: (Color online) Short time behavior of the order pa- 
rameter Am at h — ' for different system sizes A''. The 
dashed (black) line corresponds to a power law fit to the 
A'^ = 1.6 X lO'' curve. Inset: t* , defined as the time for Am 
to reach 10~^, vs. N. 
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FIG. 4: (Color on-line) A typical single-sample energy per 
spin vs. time plot, after a quench from infinite temperature 
to T = 0.99Tt for g = 24 and L = 200. Snapshots at selected 
times are also shown, with colors coding for spin values. 



Tt {[Tt - rip ^]/Tt < 10-3 for q = 5,6). It is thus very 
hard to establish whether the pseudospinodals are dif- 
ferent from the transition temperature in the thermody- 
namic limit (Binder [23| concluded that systems at least 
as large as 1000 x 1000 would be required), and the ques- 
tion of the existence of metastable phases in the thermo- 
dynamic limit has remained controversial [ill, HI] ■ 

We have determined the lower spinodal of this model 
using the short-time dynamics, as well as the pseudospin- 
odal from metaequilibrium measurements. We have done 
Monte Carlo simulations with single-spin-flip Metropolis 
dynamics on square lattices with N = L x L sites (L 
ranging between L = 200 and L = 4000) and periodic 



boundary conditions, for q = 9, 12, 24, 48, 96 and 192. 
The runs were quenches from infinite temperature, i.e. 
at constant T (typically below Tt(q)) but starting from 
a random configuration. 



A. Pseudospinodal 

For a narrow range of quench temperatures T < Tt 
the system gets stuck in a high energy paramagnetic 
metastable state, where it remains for some random time, 
after which it relaxes (relatively quickly) to the equilib- 
rium ferromagnetic state through a nucleation process. A 
typical energy vs. time curve is shown in Fig. 31 together 
with some snapshots of spin configurations, illustrating 
the nucleation process. The time for the formation of 
a critical nucleus has a log-normal distribution, whose 
average Tnuc depends both on [2^ T and L. To find 
the pseudospinodal, we look for a divergence of the re- 
laxation time of the metaequilibrium phase, so we must 
ensure that the measurements are done at times shorter 
than T^uc(T) in order to avoid entering the regime where 
domains of the stable phase have begun to grow. For this 
we consider the two-time autocorrelation function 



C(il,<2) = 



q 



1 \ N 



1 " l\ 
-J2Sia,ih),a.{t2))--j, (32) 



where ti is the time elapsed since the quench, t2 > ii, 
and the average is taken over different realizations of the 
thermal noise. Out of equilibrium, C(ti,t2) depends on 
both ti and t2 while in a stationary (meta)stable state it 
depends only on the time difference t = t2 — ti. To com- 
pute the relaxation time tr, we use correlation data only 
from the regime where it is independent of ti , thus stay- 
ing at temperatures above the metastability limit. The 
typical behavior of C (t) in this regime is shown in the in- 
set of Fig.[5]for q = 96. It is clear that the relaxation time 
is growing as one goes deeper into the metastable region. 
We defined th as the time at which C{t) falls below some 
threshold Cthr (see inset of Fig. [5]). Varying Cthr within 
a reasonable range we obtained similar results; we have 
included this arbitrariness in the error estimates. The 
behavior of tj^ as a function of temperature is plotted in 
Fig. [5] for q = 96. 

From these data we estimated the the pseudospinodal 
temperature T*. This is not an easy task: far from the 
divergence, corrections to power-law scaling can be im- 
portant, but on the other hand the pseudospinodal can- 
not be approached with equilibrium measurements. The 
estimates of T* and associated exponents will necessary 
have relatively large uncertainties. To proceed, we fit the 
relaxation time with a scaling form [iSj . 



tr^A 



T -T* 



l + B 



Tt 



(33) 



fixing c to different values. The best fit according to the 
correlation coeficcient {B? = 0.99989) was obtained for 
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FIG. 5; (Color on-line) Relaxation time vs. T/Tt for L = 1000 
and q = 96. The continuous line is a fit to the data up to 
T/Tt < 1.02 using the Eq. (|33j. Error bars are smaller than 
the symbol size. The inset shows the corresponding stationary 
correlation function vs. t for temperatures ranging from T = 
1.005 Tt (left) to T = 0.967 Tt (right). The vertical dashed 
line is T = T'. 



c — 1, for which the scahng correction to the power law 
is around 5%. However other values also give good fits. 
We find that, while the estimates of b depend on the form 
of the correction chose, the estimate of T* is within the 
interval T* = 0.95 ± 0.01. We cannot give an accurate 
estimate of b other than stating that the divergence is 
weak, with b ranging from 0.2 to 0.85. This is enough for 
our purpose of comparing T* with the spinodal obtained 
from STD, which we compute in the next section. 
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FIG. 6: (Color on-line) Short time behavior in a supercooled 
state at different temperatures for g = 96 and L — 480. Full 
lines are power law fittings for T — O.QSOTj. (a) Order pa- 
rameter m; (b) second moment m'^'. 



B. The spinodal through STD 



We now attempt to find the spinodal using the STD. 
We consider the dynamic behavior of the magnetization 
Eq. and the corresponding second moment [l^l , 

'"'^'=^i:((^g*("-.)-^)). (") 

starting from a completely disordered state. The results 
for q = 96 are shown in Fig. \E\ We observe a clear 
power law increase, m ^ and to*^^-* ^ t", spanning 
two decades in time for T = T^p^ = (0.950 ± 0.002)Tf. 
The exponents we find are rather small {0 ~ 0.06 and 
uj « 0.1). 

At variance with the mean-field case, the STD results 
show no noticeable finite-size effects. Fig.[7]shows m^^^ L? 
vs. time at the spinodal temperature Ts^^ — 0.950Ti for 
L — 480 and 960. The curves are almost indistinguish- 
able, indicating that the spinodal critical regime lasts for 
a size-independent time. Clearly, at this time correla- 
tions have not yet reached lengths of the order of the 
smallest system size. This is in contrast to the mean- 
field case (Fig. [3]), where it seems that the number of 



spins that must become correlated before faster growth 
of the stable phase begins grows with N . Again, we see 
that the span of the power law is limited, so that the val- 
ues of the exponents are not very accurate. However, our 
interest is to establish the existence of a power-law-like 
regime to define Tgp and to compare its value with the 
pseudospinodal. 

The temperature Tgp ^ identified by STD is equal to 
the pseudospinodal (here defined as the apparent diver- 
gence of the relaxation time) . The near- instability of the 

system at Tip~-* also shows up in the specific heat Ch and 
magnetic susceptibility xt- We computed these quan- 
tities in the metastable regime, where we find the same 
values either through fluctuations the energy (magneti- 
zation) or by a numerical derivative with respect to tem- 
perature (field). Plots for g = 96 are shown in Figs. [5] 
and [9l Ch and xt grow in a way compatible with a di- 
vergence at r = Tgp. However, the available equilibrium 
data span too short a range to get a reliable estimate of 
the critical exponents. 

Let us remark that that Fig. [S] is not completely equiv- 
alent to Fig. [1] (Curie- Weiss Ising). The second moment 
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FIG. 7; (Color on-line) Short time behavior of the second 
moment m'^-* at T = 0.95 Tt for different system sizes (L) 
and g = 96 Potts model. 




FIG. 8: (Color on-line) Specific heat as function of T — Tsp in 
the metastable state for q = 96 and L — 480. The full (blue) 

line is a fit with a law a(t - rip"')" ^1 + B (t - "^=p ')] ' 
The dashed (red) line is a power law fit using data for 
/Tt < 0.15. Error bars are smaller than the sym- 
bol size. The vertical dashed line is T = Ti. 



of the order parameter (part b of both figures) behaves 
the same in both cases, but not the order parameter it- 
self. This is because of the initial conditions that must 
be used. In general, at a critical point, if one starts with 
a large value of Am, one observes a power-lavir decay to- 
wards zero (the Ising spinodal corresponds to this case). 
If one starts from near zero Am, the order parameter is 
first observed to increase with a power law, then to decay 
again to zero [l^ . This initial increase is governed by an 
exponent 9, which cannot be determined from the equi- 
librium critical exponents alone (see eq. ^) ^,[121 • The 
latter is the situation in the Potts model at T = Tip"-*. In 
an equilibrium critical point, one is free to choose the ini- 
tial condition, so both situations can be observed. In the 
Potts Tsp ^ case, however, we cannot start with a high 
value of m, because that would place us automatically 
outside the metastable phase we are trying to observe. 




FIG. 9; (Color on-line) Magnetic susceptibility as func- 
tion of T — rijT' in the metastable state for q — 96 
and L = 480. The full (blue) line is a fit with a law 

a(t ~ rjp^y ^1 + b(t ~ rip"') j . The dashed (red) line 

is a power law fit using data for (t - 7;'p ') /Tt < 0.15. Er- 
ror bars are smaller than the symbol size. The vertical dashed 
line isT = Tt. 

so we can only hope to see the initial increase (the corre- 
sponding decrease at longer times is masked by the evo- 
lution to the stable phase) . In the Curie- Weiss Ising case, 
we could in principle set Am « and see something sim- 
ilar to Fig. [6^. We have failed to observe such increase, 
however, indicating that the exponent 6 is zero. Indeed, 
in equilibrium it is known that 9 = in the mean-field 
case [l^ . 

On the other hand, in the upper spinodal point of the 
the Potts model, T = Tip+^ > Tt, we can start with an 
ordered initial state, and we should again observe m de- 
crease. Indeed, we have simulated the short time dynam- 
ics for q — 96 and N = 480 taking an ordered initial state 
and tuning the temperature above the transition value 
Tt- Fig. [TUl shows the magnetization and its second mo- 
ment measured in the same way as in Fig. [T] for different 
temperatures. In Fig. llOb we can identify a nice power 
law in the second moment for Tsp'^ /Tr — 1.06±0.01 over 
more than two decades. The fitted exponent is y « 0.90. 
At this temperature. Fig. [TUk shows the magnetization 
decreases in the same way as in Fig. [1^. The main differ- 
ence is that the spinodal magnetization in the latter case 
is known exactly, while in the Potts case it is unknown 
and cannot be estimated with good precision with these 
data. Thus we do not see a power-law in the m plot, and 
we rely on m^^' to find the spinodal. 



Finally, we have repeated the above simulations for 
several values of q, computing the relaxation time and 
the STD behavior. In all cases we found that the spin- 
odal Tsp'' found by STD is compatible with the pseu- 
dospinodal T* from metaequilibrium measurements. The 
relaxation time as a function of temperature for several 
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FIG. 10: (Color on-line) Short time behavior in a superheated 
state at different temperatures for g = 96 and L = 480. Full 
line in (b) is a power-law fit for T — l.OGOTt. (a) Order 
parameter m; (b) second moment m'^-*. 



q and L = 1000 is given in Fig. [TT] (the same results were 
found for L = 2000 and 4000). This calculation was done 
using the same threshold Cthr = 0.01 for all values of q. 
We observed that for fixed T — Tgp , the relaxation time 
is non-monotonic as a function of q, with a minimum 
around g w 50 (see inset of Fig. [TT|) . and the same non- 
monotonic behavior is observed in the exponents. We 
do not have an explanation for this behavior. However, 
the growth of the relaxation time for large values of q 
is consistent with the appearance of a true singularity, 
since the g ^ oo limit is mean- field [33|. We also ob- 
served that the kinetic spinodal is strongly g-dependent 
and therefore the minimum distance to rijT'* available for 
metaequilibrium measurements increases with q. 

V In Fig. ?? we show (Tt - Tip'V^^t vs. q (we 
have included the "pseudocritical" temperatures found 
in Ref. [H for g = 5, 7). We see that rip''' systematically 
departs from Tj as q increases. In fact, the data are very 
well fitted by the logarithmic form Alog''(l -I- g — 4) with 
a = 2.81, in qualitative agreement with the behavior of 
the mean field or Curie- Weiss-Potts, and the Bethe lat- 




FIG. 11: (Color on-line) Relaxation time as function of T — 
rip^' for L — 1000 and different values of q. Error bars are 
smaller than symbol sizes. Note however that due to the error 
in the determination of Tsp, the points could be uniformly 
shifted horizontally. The continuous lines are power law fits. 
Inset: tr as a function of q for a fixed temperature T — Tap = 
0.008. 



tice solution with coordination number [3Ct l31j z = 3, as 
expected in the large q limit |32| . 



V. CONCLUSIONS 

We have shown that it is possible to define the spinodal 
point through the short-time dynamic behavior. The 
STD can be used to detect a point in the phase dia- 
gram where the dynamics is critical (albeit for a finite 
time). In mean- field systems this coincides with the 
thermodynamic spinodal defined through the vanishing 
of second derivatives of the free energy, while in finite- 
dimensional systems it serves as a definition of spinodal, a 
point where (meta)equilibrium measurements are impos- 
sible (since it is beyond the metastability limit, or kinetic 
spinodal). In the 2-d Potts model, we found that the spin- 
odal defined in this way coincides with the pseudospin- 
odal found through fitting and extrapolation of metae- 
quilibrium relaxation times. Our results are consistent 
with the scaling behavior associated with a growing cor- 
relation length. In particular, (pseudo)critical exponents 
can be measured using STD, which in the mean-field case 
we have checked with the analytical result. 

For the Potts model, this method gives a spinodal tem- 
perature different from the transition temperature at all 
g where the transition is first order, even in the thermo- 
dynamic limit. Since the spinodal provides only a bound 
for the metastability limit, this does not settle the ques- 
tion of the existence of a metastable phase, but it does 
show that the apparent convergence of spinodal and tran- 
sition temperatures is due to the extremely weak nature 
of the transition for low values of q. 

Appart from the conceptual advantage of allowing a 
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definition of spinodal points avoiding equifibration issues, 
the STD method may prove to be practically useful in es- 
tablishing bounds for metastability limits. This would be 
especially welcome for instance in systems with very slow 
dynamics such as glassforming liquids, where metaequi- 
librium measurements are out of the question (a direction 
in which some work is in progress). And even more so if 
this technique could be implemented experimentally. 

Finally, let us remark that although our results show 
that one can define a spinodal (which we call thermody- 
namic to distinguish it from the kinetic spinodal) through 
a critical-like dynamics, this behaviour is transient (ex- 
cept in mean-field), so that we do not conclude that any 
thermodynamic potential is singular at this point. The 
thermodynamic spinodal in this sense is better under- 
stood as similar to the pseudospinodal, but defined in a 
way free from the equilibration and extrapolation issues 
that pervade the pseudospinodal. In the 2-d Ising model, 
the pseudospinodal has been shown to be related to zeros 
of the partition function at complex values of tempera- 
ture and field [2l|, which approach the real plane as the 
range of the interaction is increased. We may conjecture 
a similar scenario for the Potts model, whose thermody- 
namical behavior can be determined by the zeros of its 
partition function in the complex temperature plane [20l|. 



and where evidence has been found of singular behavior 
(with divergent thermodynamical quantities) associated 
to some of those zeros [131 for g > 4. It is tempting to 
conjecture that the (pseudo)critical dynamics observed 
here is associated to such zeros, and that its lifetime is 
longer the nearer these zeros are to the real plane. This 
remains to be studied, however, but it is an interesting 
issue to investigate especially in short range interacting 
systems, where a pseudo-singularity is observed even if 
the mechanism of a vanishing free energy barrier for long 
wave length fluctuations is in principle excluded. 
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